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Abstract. We investigate the presence and properties of PAHs on the surface of circumstellar disks around Herbig 
Ae/Be stars by comparing the predictions of disk models with observations. We present results of a radiation 
transfer code for disks heated by the central star, in hydrostatic equilibrium in the vertical direction (flared 
disks). The dust is a mixture of large grains in thermal equilibrium, transiently heated small grains and PAHs. 
Special attention is given to the influence of the stellar, disk and PAHs properties on the strength of the PAH 
emission lines and their spatial distribution. The models predict an infrared SED showing PAH features at 3.3, 
6.2, 7.7, and 11.3 (jm clearly visible above the continuum, and with some of them very strong. The PAH emission, 
spatially extended, comes mostly from the outer disk region (7? ~100 AU) while the continuum emission at similar 
wavelengths, mostly due to warm large grains, is confined to the innermost disk regions (R ~few AU). We compare 
the model results to infrared observations from ISO and ground-based telescopes of some thirty Herbig Ae/Be 
stars. Most of the observed PAH features in objects with spectral type later than B9 are well described by our 
disk models and we discuss in some detail the PAH characteristics one can derive from the existing data. Objects 
with strong radiation field (generally earlier than about B9) have the 3.3 /im feature (often the only one observed) 
much weaker than predicted, and we discuss possible explanations (dissipation of the disk, photoevaporation or 
modification of the PAH properties). 
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1. Introduction 

Features of PAHs (Polycyclic Aromatic Hydrocarbons) are 
detected in a number of pre-main-sequence stars of inter- 
mediate mass with circumstellar disks, with characteris- 
tics that differ from those in the i nterstellar medium (see, 
e.g. JVan Kerckhoven et aTll2002|) . An origin of these fea- 
tures in the disk surface layers, which are directly exposed 
to the ultraviolet (UV) stellar radi ation, has be e n sug - 
gested by various authors (see, e.g.. iMeeus et all 1200 lj) . 
Recently, PAHs emi ssion has been reso l ved sp atially in HD 
10054 6, HP 97048 (Ivan Boekel et all 120041) and WL 16 



{Moore et allll998tlRessler fc Barsonvil2003(l . and shown 
to have the size typical of circumstellar disk emission. 

The presence of PAHs on disk surfaces is interesting 
for a number of reasons. First of all, as we will see in the 
following, PAHs are good tracers of the presence of very 
small particles mixed with the gas at large altitudes above 
the disk midplane, both near the star and in the very 
outer disk. Their emission can tell us if and where very 
small particles survive settling and coagulation processes 
that cause the majority of the original grain population 
to grow to very larg e sizes in the same objects (see, e.g., 
iNatta fc Tes t il l2003f) . A second point of interest is that the 



absence or presence of PAHs will have a significant impact 
on the gas physical properties and on the structure of the 
disk itself, since in photo-dissociation regions we believe 
that they contribute a large fractio n of the gas heating (via 
the photoelectric effect, s ee, e.g.. IWeingartner fc Drainel . 
l200lilHabart et alll200lh and could dominat e the H 2 for- 
mation on dust surfaces {Habart et alll2004jh 

The presence and properties of PAHs in circumstellar 
disks can be tested by comparing the observations to the 
predictions of disk models, which have been shown to re- 
produce well the g lobal spectral energy distribution (SEP) 
of the same stars {Pullemond et all 1200 lllpominik et all 
2003). We are interested in particular to check if the disk 
models can account for the observed intensity of the var- 
ious features and their spatial distribution, when known. 
To the best of our knowledge, there are no such mod- 
els available for disks around pre-main-sequence stars. 
These disks are optically thick to the stellar radiation and 
thus require a proper treatement of the radiation transfer. 
Recent models of the SEP of HP 14 1569 A (a young B9 5 
star with a complex debris disk) by iLi fc Luninel {2003) , 
which include PAHs, assume an optically thin disk, as ap- 
propriate for that particular object. 
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The paper is organized as follows. §2 provides a 
summary of the available observations of PAHs in 
intermediate-mass stars. The disk models and the radia- 
tion transfer scheme we use are described in §3.1 and 3.2. 
The dust on the disk surface is a mixture of large grains, 
which reach equilibrium with the stellar radiation field, 
very small grains and PAHs, which are transiently heated 
by the UV stellar radiation. The details of the adopted 
dust model are given in The results of the calcula- 
tions are described in §4, which also explores the effect 
of varying the stellar, disk and PAHs properties. In §5 we 
compare the model results to the observations. §6 summa- 
rizes our main results. 

2. Summary of the Observations 

From ground-based and airborne observations 
l|Brooke et all Il993t ISchutte et all Il990j) . it is known 
that aromatic infrared emission bands (AIBs in the 
following) attributed to PAHs exist in the spectra of 
some Herbig Ae/Be stars (hereafter HAeBe). About 20% 
of the HAeBe s tars surveyed to date have firmly detected 
3.3 /im feature l|Brooke et allll993|) . 

In Tabled we report for an ISO sample of HAeBe stars 
with circumstellar disks (i) the astrophysical parameters 
of each star (distance, spectral type, effective temperature, 
luminosity, and x, i-e., FUV flux at 150 AU from the star), 
(ii) the observational characteristics of the disk (evidence 
for a flared disk and disk outer radius) and (iii) the dust 
spectral characteristics (presence of AIBs, silicate bands 
and crystalline signatures). We also give, when available, 
estimate of the PAH emission extension. The strength of 
the PAH emission features at 3.3, 6.2, 7.7, 8.6, 11.2(3), 
12.7 and 16.4 fim (where most commonly observed PAH 
features occur) are given in Table Most of these sources 
are isolated HAeBe stars, with the exception of two ob- 
jects, HD 97048 and Elias 1, which are still embedded in a 
reflection nebula. These sources are intermediate between 
embedded and isolated HAeBe stars. In Tables and [3 
we also report the characteristic s of th e HAeBe stars ob- 
served at 3 /im bv lBrooke et al.l l|l993|) and of the young 
stellar object WL 16 embe dded in the p Ophiuchi c loud 
core (for a recent paper see lRessler fe Barsonvt |200"3|) . 

In the ISO sample of isolated HAeBe stars, we ca n dis- 
tinguish three different groups l|Meeus et all l200l|) : the 
first (group la) is made of stars with evidence of flared 
disk and strong or moderately aromatic and silicate bands 
(with or without crystalline signatures); the second (group 
lb) is also characterised by flared disks and presence of 
AIBs, but the silicate bands are absent; finally, in the 
third group (group II), there is no evidence of flared disk, 
weak or no AIBs, but the silicate emission features are 
strong (with various amount of crystalline silicates). The 
presence of AIBs is not correlated to any of the stellar pa- 
rameters. There is, on the other hand, an indication that 
PAH features are strong when disks are flared, since they 
are present in group I objects (both la and lb) and only 
(and rather weakly) in one source of group II. This sug- 



gests that PAH emission originates in the outer regions of 
flared disks, where dust intercepts a large fraction of the 
UV emission coming from the star. 

That in isolated HAeBe stars PAHs must be closely 
related to the star-disk system, is supported by the ob- 
servations of HD 100546 and HD 179128 (both group 
la) showing that the PAH emission occurs wi thin —150 
AU llGradv et alll200lUvan Boekel et alll2004j) and -600 
AU (|Siebenmorgen et all 120 001) from th e star, respec- 
tivelv. Moreover. Ivan Boekel et all l|2004h find that HD 
104237 (group II) is unresolved in their 10 /mr spec- 
troscopic observations using TIMMI2 at the ESO 3.6- 
m telescope with a spatial resolution of —0.9". This is 
consistent with the idea that group II sources do not 
possess a flared outer disk. As far as we know, no spa- 
tial information on the PAH emission is available for the 
other isolated HAeBe stars but their isolated nature sug- 
gests that also in these objects the PAH emission comes 
from the immediate circumstellar environment. On the 
other hand, for the two intermediate sources, HD 97048 
and Elias 1, the PAH emission probably originates in a 
significant part from the illuminating reflection nebula. 
HD 97048 observations sho w that the PAH emi s sion i s 
extended up to -1000 AU l)Siebenmorgen et all l200fj|) . 
However, recent high spatial resolution observations show 
that most of the emission of the 8.6, 11.3 and 12.7 /jm fea- 
tures cc4nes_^n_i2£t_iTOm_a_region of 200-300 AU, likely a 
disk (van Boekel et alll2004[h Additional evidence of the 
PAH emission disk origin in HAeBe stars comes from re- 
cent (sub)ar c sec-re so lution mid-IR images of W L 16 by 
iMoore et all l|l998h : iRessler fc Barsonvl l)2003h . showing 
that the PAHs emission has a size of —450 AU. Finally, 
note tha t for the HAeBe sta rs observed from the ground 
at 3 /im, iBrooke et alJ l)l993|) finds that the aromatic fea- 
tures tend to be confined to — 1000 AU. 

It is also interesting that in HD 97048 and WL 
16 the 10 /im continuum emission is spatially ex- 
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the spectral signature of silicate grains in both these 
sources (see references in Tabled, this indicates that very 
small carbonaceous grains are responsible for the contin- 
uum emission. Quite differently, the continuum emission 
in HD 100546 arises from a much smaller (<30 AU ) region 
than the PAH emission l)van Boekel et all [2004) . In this 
object, where silicate bands are present, the continuum 
may be dominated by warm large silicates grains confined 
to the innermost regions of the disk. 

Evidence of a disk origin of the PAH emission in 
HAeBe stars comes also from the extensive study of the 
ISO spectra of a sample of embedded young stellar objects, 
isolated HAeBe stars, evolved stars, reflection nebulae and 
HII regions, which reveals the presence of variations in the 
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6-9 /jm region with systematic differences between sources 
which have recently synthesized their PAHs in their ejecta 
and sources which are illumina t ing g eneral interstellar 
medium materials (jPeeters et all 12002). Combined effect 
of PAH family and anharmonicity (dependence of the peak 
position with the temperature of the emitters) could be at 
the origin of these variations. The 6.2 and 7.7 /zm fea- 
tures in isolated HAeBe stars which resemble those in 
evolved stars are very different from those in sources with 
ISM material. Moreover, large variations in the ratio of 
the band strengths in the 10-15 /im r egion are observed , 
especially in the 11.2/12.7 [im ratio ijHonv et all l200l|) . 
Together with the good correlation found between the 11.2 
and 3.3 um band and between the 12.7 and 6.2 /xm band, 
IHonv et all l|200ll ) suggest that the spectra of HII regions 
where the 11.2 and 12.7 /tm are typically equally strong 
are due to small or irregular (and preferentially ionized) 
PAHs. On the other hand, the spectra of evolved stars with 
a small 12.7/11.2 /an ratio may arise from large (~100-150 
C-atom) compact (and preferentially neutral) PAHs. The 
12.7 /im band cannot be seen in the low resolution spec- 
tra of isolated HAeBe stars and the 12.7/11.2 /mi ratio 
is not know n. But recent ground-based o bservations (see 
reference in IVan Kerckhoven et all 12002^ confirm that, if 
present, the 12.7 band must be very weak compared 
to the 11.2 fim. This indicates that the PAHs are more 
large and less ionized in isolated HAeBe sources. 

In summary, there are indications that PAHs in HAeBe 
stars are somewhat different from those in the ISM, and 
more similar to PAHs in evolved stars. However, it is dif- 
ficult to fully characterize the PAHs properties in HAeBe 
stars from what is known at present, and in the following 
we will start our analysis from the assumption that PAHs 
in the ISM represent a good first approximation to the 
properties of PAHs in disks as well. We will then discuss 
possible differences, and come back to this point in our 
conclusions. 

3. Models 

In this section, we describe the disk models and the ra- 
diation transfer used and give the details of the adopted 
dust model. 

3.1. Disk Models 

We consider disks heated by irradiation from the central 
star, in hydrostatic equilibrium in the vertical direction, 
with gas and dust well mixed (flared disks). We com- 
pute the disk structure us ing an improved version of the 
Chiang fc Ooldr eicbl lll997h two-layer models, described in 
Dull emond et alJ 1)20011) . For a given star, the disk struc- 
ture (i.e., pressure scale height H p and flaring angle a) is 
completely defined once the inner and outer radii (Rj n and 
Rd), the surface density distribution £ = £o(i?/i?o)~ p 
with i?o a fiducial radius, and the dust model are speci- 
fied. 



In these models, the stellar flux impinging with the 
flaring angle a into the disk is absorbed in the upper layers 
of the disk, which will reradiate half of the flux away from 
the disk and half down into the disk's deeper layers. Here, 
we compute the structure of the disk self-consistently us- 
ing an iterative scheme that gives the flaring angle a and 
the midplane temperature at each radius. The radiation 
transfer is solved in an effective, simplified way, which, 
however, gives a disk structure very similar to that ob- 
tained with a full, much more time consuming radiation 
transfer treatement l|Dullemond fc NattaLl2003|) . 

The models are appropriate for disks that are optically 
thick to the stellar radiation. This is generally the case for 
disks around pre-main-sequence stars, up to very large 
radii (e.g., 5000 AU for a disk mass M d ~ 0.2 M , p = 1). 

3.2. Radiative Transfer 

In spite of the success of the two-layer codes in predict- 
ing the disk structure and the overall SED, they tend 
to overestimate the disk e mission in the near-infrared 
<|Dullemond fc NattaL [2003). We have therefore decided 
to implement a full ID radiation transfer scheme to com- 
pute more accurately the emission in the wavelength range 
where the PAH features occurr. In practice, we first com- 
pute the disk structure using our two-layer code, with the 
adopted dust model (which includes transiently heated 
species), then, keeping fixed the disk structure, we com- 
pute the emerging spectrum using the radiation transfer 
code. This procedure is not fully self-consistent, but we 
have checked that it does not introduce any significant 
inaccuracy in our results. 

Dust scattering is incorporated. Scattering of stellar 
light may be relevant because it can lead to a substantial 
reduction of the flux available for heating the disk and be- 
cause it allows stellar photons, after one scattering event, 
to penetrate much deeper into the disk then one would 
expect given the smallness of the flaring angle a. 

To calculate the emission of transiently heated small 
grains and of PAHs, we com pute the temperature distri- 
bution functions dP/dT as in ISiebenmorgen et alJ il 992) 
follow ing the method outlined bv lGuhathakuTtXfcpTanie) 
Jl989l) . 

To determine the emission from the disk, we divide it 
into rings of radius r and small radial width Ar. At the 
center, at r — 0, sits the star with parameters M+, L+, R* 
and T e ff. For each ring, we compute, as described in the 
appendix, the radiative transfer perpendicular to the disk, 
in z-direction, and then add up the contributions from all 
rings. 

3.3. Dust model and PAH properties 

The properties of dust in circumstellar disks differ signif- 
icantly from those in the ISM (see and from star to 
star. It is difficult to build a general, realistic dust model. 
However, this is not the purpose of this study, where we 
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only want to test the hypothesis that the observed PAH 
features have their origin in the surface layers of circum- 
stellar disks. Therefore, we will start with a dust model 
that includes large grains (i.e., grains large enough to be 
in thermal equilibrium with the local radiation field) with 
properties that roughly account for the observed SEDs (in- 
cluding the silicate emission at 10 /zm) of most pre-main- 
sequence stars. To these, we will add very small particles 
and PAHs, with abundances and properties typical of the 
ISM, that we will then vary. 

a) Large grains. The large or big grains (BGs in the follow- 
ing) consist of grap hite and silicates with optical constants 
from lDraind 1)19851) . These grains have a MRN size distri- 
bution (n(a) oc a -3 ' 5 ) between a minimum and a maxi- 
mum radius ai =0.01 /im and a u =0.74 /jm for silicates 
and ai =0.01 /im and a u =0.36 /im for graphite. We as- 
sume that all silicon and about 2 /3 of the cosmic carbon 1 
([C / H]bgs = 1.5x 10 -4 ) are in these grains. This prescrip- 
tion gives a good fit to the mid-IR SED (whi ch originates 
from the disk surface) of many HAe stars l|Natta et all 
200lJ; we have found in our modeling that the details are 
not very important for the purpose of this paper. 

b) Very small grains. To the BGs population of grains, we 
add a tail of very small graphite particles (VSGs in the 
following) with the same size distribution n(a) oc a -3 ' 5 
betw een a; =10 A to a u =100 A. Optical constants are 
from lDrainel l(l985h . The fraction of carbon in VSGs is 15% 
that in BGs ([C/H] VSGs = 0.22 x 10~ 4 ). VSGs and BGs 
graphite grains probably form a continous population. The 
division between BGs and VSGs is maintained mainly for 
computational purposes. Note t hat there is no evide nce in 
the ISM of very small silicates l|Li fc Drainel l2001bl) . 



c) PAHs. Finally, we include in our dust model a popula- 
tion of PAHs. In the ISM, PAHs are made up of a few tens 
up to a few hundreds carbon atoms; for reasons of simplic- 
ity, we start our models with one PAH size, Nc = 100. 
The hydrogen to carbon ratio is H/ C = f H x (6/N g)° 5 
(case of compact symetric PAHs, see lOmontl Il986|) with 
fn the hydrogenation fraction of the molecule. We start 
our models assuming essentially fully hydrogenated PAHs, 
i.e, fjj = 1- The carbon locked up in PAHs has an 
abundance of [C / H]pah — 5 x 10~ 5 inferred from the 
12 / im emission per hyd r ogen in typical Galactic cir- 
rus (iBoulanger fc PeraultL Il988j) and photo-dissociation 
regions l|Habart et all l200l[) and from comparison be- 
tween observations of dust galactic e mission and extinc- 
tion with detailed model calculations llDesert et all FlQQO: 
iDwek et all Il997t iLi fc Drainel EoOlal). We take the ab- 
sorption cross section as defined bv lLi fc Drainel l|2001a|) 
based both on laboratory data and astrophysical spectra. 
They consider the bands at 3.3, 8.6, 11.3, 11.9 and 12.7 



/im from vibrational modes of the aromatic C-H bond; the 
strong bands at 6.2 and 7.7 /*m due to vibrations of the 
aromatic C-C bonds; and a few weak features probably 
caused by C-C bending modes at 16.4, 18.3, 21.2 and 23.1 
/im. The band profiles have a Drude shape; the integrated 
cross sections band profiles, positions and widths of the 
lines are found in their Table 1 2 . The cross section in the 
visible-UV is determined from the laboratory absorption 
spectrum of small species with a cut-off in the visible to 
ne ar-IR whose wavele ngth increases with Nc, as defined 
bv lDesert et"aH l|l990|) . 

With respect to the charge, we start assuming that 
PAHs are mostly neutral. Ideally, one can calculate the 
steady state charge distribution of PAH from the bal- 
ance between electron capture rates and the photoelec- 
tron emission rates plus the ion capture rates (see, e.g., 
IWeingartner fc Drainel l200l|) . Since little is known re- 
garding the electron density n e and its distribtion in the 
surface layers of disks, precise determination of the PAH 
charge distribution in disks is not possible. Instead, we 
can estimate the ratio between the photoionisation rate 
of the grain to the electron-grain recombination rate, 

7, which can be approximated by 70 x f(Nc) N^/ 3 x 
1/2 

X Tgls /n e with f(Nc) a function giving the dependence of 
the photo-electron ejec tion probability with the size (see 
Bak es fc TielensL Il994f) . For a disk heated by a typical 
HAe star, we have \ ~ 10 5 an d n-H ~ 10 7 cm~ 3 on the 
surface disk layer at 150 AU from the star (see next sec- 
tion). Thus, assuming n e /nn ~ 1.5 10~ 4 (as in typical 
photo-dissociation regions of the ISM) and T gas ~ 100 K, 

we find x Tgls /n e ~ 500 and for a PAHs with N c ~ 100 
we estimate 3 that 7 is inferior than 1. This implies that 
PAHs must be neutral or negatively charged in the outer 
disk regions. In the inner disk region and for stars with 
higher x, however, 7 becomes higher and PAHs will be 
in part positively charged. To keep the model simple, we 
neglect that we probably have a mixture of neutral and 
charged PAHs in disks and we assume that PAHs are neu- 
tral. The effect of charge is briefly discussed in M4.2.1I 

Upon absorption of an energetic photon, small PAHs, 
with an insufficient number of internal vibrational modes 
in which to distribute this photon energy may dissociate. 
In a strong FUV radiation field (x > 10 2 ), destruction 
occurs during multiphoton events where PAHs absorb an 
energy of more than 21 eV in a time shorter than the cool- 
ing t ime ( Guhathakurta fc Drainel Jl989t ISiebenmorgenl 
Il993h . Recently. ILi fc Limine! ( |2003h estimated the pho- 
todestruction rate for the HD 141569A disk and suggest 
that PAHs with N c < 40 - 50 will be photolytically unsta- 



1 ISnow fc Witt! lll995l) have discussed the carbon abundance 
in stars newly formed from the ISM and conclude that the 
carbon abundance in the local ISM is (2.24±0.5) x 10~ 4 . 



2 In order t o fit th e observed IR spectrum of the diffuse ISM, 
ILi fc Drainel feOOlaf) include enhancement factors Ee.2, -B7.7 
and -Eg. 6 in the integrated cross sections for the 6.2, 7.7 and 
8.6 fim bands. For the template model, we adopt Eq,2, -E7.7, 
E&.e equal to 1, which recover typical laboratory values, but we 
will also consider model with Ee.2, -E7.7, i?s.6>l (see H4.2.51 . 

3 To estimate 7, we tak e 70 ~ 3 X 10 -6 and f(Na) ~ 11 as 
in lBakes fc Tielensl Jl994l) . 
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Fig. 1. Absorption coefficient (in cm 2 per gram of gas) of 
large grains (BGs; dotted line), very small grains (VSGs; 
dashed line) and PAHs of 100 C atoms (solid line). See 
text for details. 

ble in the inner R < 100 AU region around the star (where 
X > 10 5 ) during the life time of the disk. For a PAH with 
Nc ~ 40 — 50, the maximum temperature reached just 
after the deposit of an e nergy of 21 eV is about 2000 K 
(see, e.g., lOmontl Il986l) . In our model, we assume that 
PAHs are destroyed when the probability that they are at 
a temperature T > T evap exceeds a critical value p e vap, 
i.e., when P(T)dT > p eV a P - Such procedure is nu- 

merically simple and can probably correctly evaluate evap- 
oration. However, we do not derive p eva p from thermo- 
dynamical considerations. We use T evap = 2000 K and 
Pevap — 10~ 8 and checked that for this choice, in the case 
of small PAHs, the results are similar to the assumption 
that the grains are photo-destructed if they absorb an en- 
ergy of more than 21 eV in a time interval shorter than 
the cooling time. 

In summary, we adopt a dust model with BGs, VSGs 
and PAHs. The silicate abundance in dust is [Si/H] = 
3 x 10~ 5 , and the total carbon abundance in dust is 
[C/H] = 2.22 x 10~ 4 . Of this, 23% is in PAHs, 10% in 
VSGs and 67% in large grains. Fig. ^ shows the adopted 
absorption coefficient for the various species. For a radi- 
ation field that has a blackbody distribution at T=10 
K, BGs will absorb 50% of the radiation, VSG 11% and 
PAHs 39%, respectively. The properties of BGs are some- 



what different from those of grains in the ISM, and have 
been chosen to approximately fit the properties of dust 
on disk surfaces. VSGs and PAHs have properties derived 
from the ISM. In the following, however, we will explore 
the effects on the predicted spectra of some of our assump- 
tions for the PAHs properties. Note, however, that we will 
not try to reproduce the peculiarities of the features ob- 
served in some isolated HAe stars, such as the shift of 
the peak positions, but we will limit our discussion to the 
band strengths. 

4. Results 

In this section, we describe the results of the calculations 
and explore the effect of various stellar, disk and PAHs 
properties. 

4.1. Template Model 

Firstly, we discuss the results for a model that we will use 
in the following as a template. 

The disk is heated by a typical HAe star with effective 
temperature T e ff=10 500 K, luminosity L*=32.0 L Q and 
mass M*=2.4 M Q . At a distance of 150 AU, the FUV 
radiation field integrated between 912 and 2050 A is x = 
10 5 time s the average in terstellar radiation field (1.6 x 10 -6 
W m- 2 : lHabimlll96S|) . 

The disk mass is ~ 0.1 M©. Its inner radius is Ri ra = 
0.3 AU, close to the BG evaporation radius, the outer 
radius is Rd=300 AU. The surface density profile is £ = 
S x (r/Ri)^ 1 with £ = 2 10 3 g cm -2 . We assume that 
only half of the star is visible from any point on the disk 
surface = 0.5). Fig. |3 shows the run with radius of the 
pressure (H p ) and photospheric (H s ) scale height and of 
the flaring angle a (bottom panel). In the top panel we 
show the temperature of silicates and graphite grains of 
size a = 0.01 /im at H s . The gas density at H s decreases 
from 7 x 10 9 cm" 3 at R m to 3 x 10 6 cm" 3 at R d . 

Each side of the disk intercepts and reprocesses about 
10% of the stellar luminosity. The fractional contribution 
to the reprocessed radiation of BGs, VSGs and PAHs are, 
respectively, 60%, 13% and 27%. Note that the fraction 
emitted by PAHs is lower than the fraction of the stellar 
radiation they absorb (about 39%, see precedent section) 
because of their evaporation in the inner part of the disk. 
About 10% of the total absorbed flux is lost in PAH evap- 
oration. 

Fig. El shows the calculated spectrum in the 2-40 /jm 
range of a star/disk system seen face-on at a distance of 
150 pc. PAH features, which form in the optically thin 
disk surface, are clearly visible, and some of them are very 
strong; in Table we give the integrated flux in the 3.3, 

6.2, 7.7, 8.6, 11.3, 12.7 and 16.2 features after contin- 
uum subtraction. 

The continuum is mostly due to the large grains which 
are very hot in the inner regions (see Fig. |2J). It is possi- 
ble to see, for example, the broad feature due to silicate 
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Fig. 2. Disk properties for the template model. The lower 
panel shows the run with radius of the pressure scale 
height Hp, the photospheric scale height H s and the flare 
angle a divided by 20 deg for convenience, a is 1.4 deg 
at the inner radius, and reached 7 deg at the outer disk 
radius. The top panel shows the temperature for graphite 
(solid line) and silicate (dashed line) grains of size 0.01 
/zm at H s . 



Fig. 3. Spectral energy distribution of a star/disk system 
(template model) seen face- on at a distance of 150 pc. The 
solid line shows the emission of the disk, the dotted line 
that of the stellar photosphere. The dashed line shows the 
SED of a model which has all the parameters of the tem- 
plate model but no VSGs and no PAHs. The stellar param- 
eters are T cfi = 10,500 K, L*=32 L Q , M*=24 M Q . The 
disk is flared, it has a mass of 0.1 Mq, an inner radius of 
0.3 AU and an outer radius of 300 AU (see Fig. QJ). 



emission peaking at about 10 /zm under the much nar- 
rower PAH features. Some contribution to the continuum 
in the wavelength range A < 10 /zm is also due to the 
VSG, as shown by the SED of a model which has all the 
parameters of the template model but no VSGs nor PAHs 
(dashed line in Fig. ^ . 

The PAH emission is much more extended than the 
adjacent continuum, as shown by two figures, Figs. 0] and 
The Fig.0]plots as function of the projected radius the 
cumulative intensity of the feature and continuum at 3.3, 
6.2 and 11.3 fim. The continuum reaches 50% of its in- 
tensity at a very small radius (typically 2-5 AU) while the 
feature does so at larger radii (about 30 AU for the 3.3 /zm 
feature and 80 AU for the 6.2 and 11.3 /zm features). This 
behaviour basically reflects the different excitation mech- 
anism of PAHs, i.e., single-photon excitation, and of the 
continuum at the same wavelength, which is dominated 
by the emission of BGs, in thermal equilibrium with the 
radiation field. Only very hot BGs, located in the inner- 
most disk, emit at 3 fj,m whereas PAHs also emit at the 
outer edge of the disk in the 3.3 /jm feature. 

The ratio of short wavelength over long wavelength 
band strength decreases with the intensity of the FUV ra- 
diation also for transiently heated particles. Fig. 01 shows 
that the contribution of the outer disk to the 3.3 /zm fea- 



ture is small, while it is much larger for the features at 
6.2 and 11.3 /urn. We plot for comparison the run with 
radius of the cumulative fraction of the FUV radiation in- 
tercepted by the disk surface. The similar behaviour with 
the cumulative intensity of the less energetic PAH features 
reflects the fact that the emission of PAHs roughly scales 
with the intensity of the FUV radiation field (see H4.3.1I 
for more details). 

Fig. shows the intensity profile obtained by convolv- 
ing the computed intensity with a two-dimensional gaus- 
sian of FWHM=0.1" . As expected from Fig.0J the features 
at short wavelengths are strong in the inner hot part of 
the disk decreasing rapidly in the outer cold part, whereas 
those at long wavelengths are weaker in the inner region 
(mostly because the continuum is higher) but more ex- 
tended (the temperature required to excite them is lower). 
In other words, the 3.3 /zm probes the inner (R < 100 AU) 
region of the disk while the other features probe its outer 
parts. 

4.2. PAH properties 

In this section, we investigate the effect of the charge, 
hydrogenation parameter, size, photoevaporation and ab- 
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Fig. 4. Cumulative intensity of the feature ( continuum 
subtracted, solid lines) and continuum (dotted lines) at 
3.3, 6.2 and 11.3 \im as function of the projected radius. 
We plot for comparison the run with radius of the cumula- 
tive fraction of the FUV radiation intercepted by the disk 
surface (dashed lines). 



Fig. 5. Surface brightness profiles convolved with a two- 
dimensional gaussian of FWHM—0.1" at the peak of the 
feature (no continuum subtracted, solid lines) and in the 
adjacent continuum (dotted lines) for the 3.3, 6.2 and 11.3 
/im features. In the top panel, we also show the photo- 
spheric stellar emission at 3.3 \im convolved with the same 
beam (dashed line). As in Fig.\Q the star/disk is seen face- 
on at a distance of 150 pc. 



sorption cross section of PAHs on the predicted intensity 
of the features. 

4.2.1. Charge state 

In Fig. El we show the spectra and the brightness distri- 
bution for the same model shown in Figs. Eland but for 
PAH cations. The PAH spectrum is a strong function of 
the PAH charge: the strength of the C-C stretching modes 
between 6.2 and 8.6 /iin is much stronger for ionized PAHs; 
the 7.7 /im C-C feature to the 11.3 /im C-H feature ra- 
tio is increased by a factor of ^5 (see Table OJ. Further, 
whereas neutral PAHs show a strong 3.3 /im C-H features, 
the band almost vanishes when the PAHs are positively 
ionized. 

The IR properties of PAH anions closely resemble 
those of PAH cations except for the very strong 3.3 
um enhancem e nt in t he anion llSzczena nski et all Il99fit 
iHudgins et all HoOO). However. iBauschlicher fc Bakes! 
(2000) predict that PAH anions have band strengths in- 
termediate between those of neutrals (strong 3.3 /im C-H 
stretching and 11.3, 11.9 and 12.7 /im out-of-plane C-H 
bending modes) and PAH cations (strong 6.2 and 7.7 /im 
C-C stretching and 8.6 /im C-H in-plane bending mode). 



4.2.2. Hydrogenation parameter 

Another effect which can significantly affect the PAH spec- 
tra is the dehydrogenation. Strong dehydrogenation would 
make the C-H features disappear. However, theoretical 
studies of the dehydrogenation in the ISM indicate that 
dehydrogenatio n is relevant only for the smalles t PAH s 
with N c <25 llTielens et all Il987t lAlla.in et all llflflrf) . 
Upon absorption of an energetic photon, large PAHs will 
not have C-H bond rupture since the absorbed energy 
will promptly be redistributed among many vibrational 
modes. Nevertheless, in the inner disk regions and for stars 
with high UV fields dehydrogenation could become impor- 
tant even for large PAHs. To estimate this effect, we run 
a model with PAHs partially dehydrogenated, i.e., hydro- 
genation fraction /#=0.5. The results are shown in Fig.[HJ 
The C-H features (3.3, 8.6, 11.3, 12.7 /im) become lower 
(by a factor ~2, see Table EJ, and the C-C features (6.2, 
7.7 /im) become slightly higher (as expected from energy 
conservation) . 

4.2.3. Size 

We now look at the influence of the PAH size. In Fig.0 we 
show a model for smaller PAHs with Nc=4:0. The relative 
strength of the PAH emission bands depends strongly on 
the PAH size : small PAHs radiate strongly at 3.3 /im, 
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Fig. 6. The top panel shows the SEDs of models contain- 
ing positively ionized PAH (solid thick line) and partially 
dehydrogenated PAHs with fn = 0.5 (dotted line). All 
the other parameters are as in the template model. For 
comparison, we show again also the SED of the template 
model (solid thin line). Note that the fluxes of the ionized 
PAH model and of the template model have been mul- 
tiplied by factors 10 and 100, respectively. The bottom 
panel plots the brightness profiles convolved as in Fig. 
at the peak of the 3.3, 6.2 and 11.3 /im features (no contin- 
uum subtracted) , as labelled. The lines refer to the ionized 
PAH model (solid thick) , the dehydrogenated PAH model 
(dotted) and to the template model (solid thin). 
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Fig. 7. Same as Fig.\Qfor model containing neutral PAHs 
with Nc=40 (solid thick lines) and ionized PAHs with 
Nq=40 (dotted lines). Again, we show also the template 
model by solid thin lines. 

while larger (colder) PAHs emit most of their power at 
longer wavelengths. The 3.3/6.2 /im ratio is ~ 4 times 
higher for small PAHs than for large PAHs (see Table |3J). 
In addition, since the H/C ratio decreases with the size 
(see H3.3[l . the C-H features increases for small PAHs. 

We note that decreasing the PAH size goes generally 
in the opposite direction than increasing the fraction of 
ionized PAHs: the emission in the 6-9 /jm range decreases 
for smaller PAHs while for ionized PAHs it increases (see 
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H4.2.1I) . On the other hand, the emission at 3.3 and 11.3 
/im goes up when PAHs become smaller, whereas it goes 
down when the fraction of ionized PAHs increases. 

Therefore, we expect to find for small and ionized 
PAHs results similar to those predicted by the template 
model with large and neutral PAHs. For illustration, we 
show in Fig. [3 a model with small (Nc = 40) and ionized 
PAHs. In this case, we find in fact that the strength of the 
6.2, 7.7 and 11.3 /im features are roughly similar to those 
predicted by the template model (see Table |2J. However, 
size and ionization do not compensate each other fully 
and the 3.3 /im band strength is considerably lower (by a 
factor of 10) than in the template model. But note that 
for small PAH anions, as mentioned in H4.2.1I the 3.3 /im 
feature strength will be higher. 

Finally, we find that for small PAHs which are mostly 
destroyed in the inner (hot) part of the disk - R < 70 AU 
for iVc=40 - the brightness profiles of the features are af- 
fected. The emission is decreased in the inner part and a 
plateau appears reflecting the compensation between the 
graduated coming out of PAHs and the decrease of the 
FUV flux. In the next section, we discuss the photoevepo- 
ration effect. 

4.2.4. Evaporation 

With the formalism adopted here, we find that large PAHs 
with Nc = 100 will not be photo-evaporated in most of 
the disk, except in the R > 20 AU region, while smaller 
PAHs with A r c=40 will be mostly photo-evaporated in 
the R < 70 AU region (as mentionned above). Larger 
PAHs are more stable since they can easily accomodate 
the absorbed photon energy. But PAHs evaporation, as 
discussed in H3.3I is a complex process, and it is possi- 
ble that our treatment of it underestimates its effects. For 
stronger photoevaporation, the strength of the most ener- 
getic features which comes mostly from the inner region 
will decrease and their brightness distribution will be af- 
fected; the effect on the bands at longer A coming from 
the outer region will be on the contrary negligible. 

On the other hand, it is possible that continuous 
replenishment of PAHs via sublimation of icy mantles 
(in which interstellar PAHs may have condensed during 
the dense molecular cloud phase) or by acc reting carbon 
atom s and/or ions from the gas (see, e.g.. lAllain et all 
Il996h could maintain PAHs throughout all the disk. For 
the template model without evaporation, the strength of 
the 3.3 /im feature will increase by a factor of ^2 while 
the other features at longer wavelength will not. 

4.2.5. Absorption cross section 

Finally, we discuss the influence of the absorption cross 
section. First, the PAH spectra depend on the integrated 
cross sections of individual features, cr}^', which, unfortu- 
nately, differ among various lab measurements or between 
experimental and theoretical studies. To get an idea of how 



these uncertainties affect the PAH spectra we have com- 
puted models with enhancement factors of 2-3 in cr}^. In 
particular, we have taken for the 6.2, 7.7 and 8.6 /im bands 
E R 2=3 , _E 7 . 7 =2 and E s ^ 6 —2 as suggested bv lLi fe Drainel 
l|2001aj) in order to fit the 6.2/7.7, 7.7/11.3 and 8.6/7.7 /im 
band strength ratios observed in ISM spectra. Compared 
to the template model, the 6.2/7.7 and 7.7/11.3 /im ratios 
are higher by a factor ~2. 

Second, the intensities of the PAH bands depend on 
the absorption cross section in the visible and UV which 
are rather uncertain. The cross section (cr^y) defined in 
t |3.3l has been derived for planar PAHs fro m laboratory 
measurements of small species (N c < 30. IJoblin et all 
Il992|) an d, for bigger species , from optical constants for 
graphite l|Verstraete fc Leeerl [l992). However, for astro- 
nomical PAHs with more complex 3-dimensional struc- 
ture very little is known. To estimate the effect of this un- 
certainty we have com puted a model with a ^R r for small 
graphite spheres as in iDraine fc Lei l|l984h . The power 
absorbed by Ojjy is ^2 times lower than for a^Jy and, 
consequently the intensity of the features is reduced by 
-2. 

4.3. Star and Disk Parameters 

The SED of the template model depends not only on the 
dust model, and in particular on the PAHs properties, but 
also on the star and disk properties. In this section, we 
investigate the influence of star/disk system parameters 
on the PAH emission features. 

4.3.1. Stellar radiation field 

The range of values of the FUV radiation field is large 
even within the rather limited range of spectral types we 
investigate in this paper, with \ (the FUV field at a dis- 
tance of 150 AU from the star) varying from 10 6 to few 
10 2 (see Tabic \ changes because both the luminosity 
and the effective temperature of the stars change. We have 
therefore computed a set of disk models for stars of dif- 
ferent spectral types assuming that they are on the main 
sequence, and plotted the results in Fig.|Hlas function of 
the corresponding value of x- The spectral types are also 
displayed in the figure. 

As expected, the strength of the PAHs features in- 
creases with x- The dependence is not linear (see Fig. 
|8J , because the fraction of the stellar radiation absorbed 
by PAHs depends on its wavelength dependence, i.e., on 
T c g. For example, the power absorbed by PAHs, nor- 
malised by x, is t wo times higher for T e ff=10 000 K than 
for T cff =15 000 and ten times lower than for T off =6 000 K. 

Fig.|Hlshows the intensity of the 3.3, 6.2 and 11.3 /im 
features as function of x- One sees that the 3.3 /im one in- 
creases faster than the others, since it is the most sensitive 
to the PAH excitation temperature. One can also note that 
the competing effect, i.e., the fact that PAH evaporation 
is larger for higher x, reducing the 3.3 /im feature more 
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Fig. 8. Strength of the PAH features at 3.3 (solid line), 6.2 
(dotted line) and 11.3 /im (dashed line) for ZAMS stars. 
The scale on the upper axys shows the spectral type, that 
on the lower axys the corresponding value of \. The thin 
solid line shows a linear dependence Ipah oc %. Disk and 
dust parameters are as in the template model, which has 
X = 10 5 (L*=32 L Q and T cff =10 500 K). 

than the others, does not compensate entirely the effect of 
higher excitation temperatures. However, one should keep 
in mind that for smaller PAHs (Nc < 100), the effect of 
evaporation could become much larger. 

4.3.2. Disk geometry 

The dependence of the PAH emission on disk parameters 
is easily understood. The quantity that most affects the 
resulting spectrum is the disk flaring angle, which deter- 
mines at each radius the fraction of FUV intercepted by 
the disk surface. Lower values of the flaring can be caused 
by a variety of reasons, for example if the dust settles to- 
ward the disk midplane. Also, if the disk mass is fixed, a 
steeper dependence of the surface density profile than p=l 
tends to make the outer disk less flared, as the ratio of the 
photospheric to the pressure scale height increases. In all 
cases, the value of the flaring angle in the inner disk is 
pract ically fixed, since it is dominated by its geometrica l 
part ^Chiang k Goldreichl Il997t iDullemond et alll200l|) . 
Therefore, less flared disks have lower emission in all fea- 
tures, but a higher ratio of the 3.3 /jm intensity, which 
forms mostly near the star, where the flaring is domi- 



nated by the geometrical part, to that of the 6.2 and 11.3 
/jm ones, which have a much larger contribution from the 
outer disk. In the extreme case of a geometrically flat disk, 
the PAH feature strenghts decrease by 3 order of magni- 
tude but the ratio 3.3/11.3 /im increases by a factor of 3 
and only the 3.3 /(m band not completly disappears. The 
3.3 /im band to continuum ratio decreases by a factor of 
4 and becomes equal to 1.1. 

In a similar way, when the disk is larger, all features get 
stronger, but the 3.3 /urn band less than the others. If FLi 
becomes very large, the disk turns optically thin to the 
FUV radiation and the PAH emission does not increase 
further. 

Another parameter to be considered is the disk incli- 
nation with respect to the observer. If the disk is inclined 
it will emit a lower continuum emission at long wave- 
length; for an inclination of 9 = 60 deg (9 = for face-on 
disks) the models predict a continuum at A >10 /im about 
20% lower than face-on models, while at shorter wave- 
lengths there are practically no changes. The continuum- 
substracted intensity of all the PAH features we compute 
will not be affected at all. 

A major source of uncertainty in our results derives 
from the lack of knowledge of the detailed structure of the 
inner disk. We have assumed in our calculations that the 
disk is truncated at the dust sublimation radius, and that 
each disk face is illuminated by 1/2 of the stellar surface. 
In fact, there is good evidence that in most HAe stars 
the disk has a puffed-up rim at the dust sublimation ra- 
dius, due to the direct illumination of t he ste llar radiation 
l|Natta et all l200li IDullemond et all l200lh . The struc- 
ture of this rim, and the effects of the shade it projects 
further out on the disk are rather complex, depending on 
the detailed 2D radiation transfer for a mixture of grains 
of different properties. The correct treatment of the ef- 
fects of this region on the PAH emission is well beyond 
the current capabilities of our models. However, we can 
make some qualitative estimates in the following way. A 
puffed-up inner rim will emit a strong continuum at short 
wavelengths. This affects mainly the ratio of the 3.3 /im 
band to the continuum; considering that models that in- 
clude the rim predict a 3.3 tt m continuum a fact or of ~ 3 
larger than our models (see iNatta et al. I l200ll> . the 3.3 
/xm peak to continuum ratio will decrease by a factor of 
2. However, the continuum-subtracted intensity of the 3.3 
/tin feature will not be affected; also, the effect on the 
other features at longer wavelength is negligible. 

The geometrical effect of the rim on the disk illumina- 
tion can be much smaller than we have assumed by taking 
0* = 0.5, to the limit of being entirely negligible. In this 
case, each point on the disk will be illuminated by the 
entire stellar surface, and the fraction of intercepted stel- 
lar radiation will increase by a factor of 2 with respect to 
our template model. We have checked that this results, 
with a large degree of accuracy, in a factor of 2 increase 
of all the strong PAH features. The total IR excess will 
also be a factor of 2 larger, going from about 20%, for 
a face-on disk, to about 40%. Most HAe stars have in- 
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deed IR excesses within th is interval l|Meeus et all l200ll 
Dull emond fc NattaL |2003) , and one expects that more re- 
alistic models will be roughly intermediate between these 
two cases. However, it should be noted that one can expect 
also some modest change in the ratio of features at differ- 
ent wavelengths, as the rim shade will affect differently 
different regions of the disk. 

5. Comparison with the observations 

A first test of our models can be made by comparing the 
predicted intensity of the most commonly observed PAH 
features (i.e., at 3.3, 6.2 and 11.3 /jm) to the observa- 
tions. We do that in Fig. [§| which plots for each feature 
the intensity normalized to the FUV flux of the star, as 
function of the value of x (i- e -> the FUV field at 150 AU 
from the star). Both quantities are distance- independent. 
The solid lines show the prediction of models for main- 
sequence stars of different spectral type (and x)> as de- 
scribed in §4.3.1. The lower line refers to models where 
0* — 0.5, the upper lines to models with (p* = 1. The thin 
lines show an uncertainty strip of ± a factor of 2 around 
these models. This, as we have discussed, is a conservative 
estimate of the uncertainties due to deviations from the 
assumptions of the template models, in the disk descrip- 
tion (different surface density profiles, flaring etc.) and in 
the assumed dust model. 

5.1. PAHs in disks: consistent with observations? 

In the following, in order to see if PAHs in disks are con- 
sistent with observations, we check for each star group 
described in Sect. [21if the disk models can account for the 
observed PAH intensities. 

5.1.1. Flared disks 

O bjects with flared d isks, according to the classification 
of iMeeus et all (|200l[) ( group I) are marked as crosses in 
Fig. This group dominates among detections, as ex- 
pected from our models. They are in general well ac- 
counted for by disk models with "standard" dust, as 
defined in H3.3I This is true also for the two younger 
sources, i.e., HD 97048 and Elias 1 (shown by trian- 
gles), especially if we correct for contamination in the two 
longer wavelength bands by the associated reflection neb - 
ula (about ~40% for HD 97048: Ivan Boekel et alll2004h . 
We also show by as terisks the HAeBe stars observed by 
IBrooke et al ] d!993h at 3.3 /tm. They cover a range of x 
from about 10 4 to very high values (~ 10 6 ), and we do 
not have detailed informations on their disk properties. 
However, the objects with x ^ 10 4 are consistent with our 
disk predictions. No upper limit in objects with x ^ 10 5 
is significantly lower than the model predictions. 

In 8 of these objects one has a measurement (or signif- 
icant upper limit) of both the 3.3 and the 6.2 /mi features. 
The ratio ranges from 0.1 to about 0.5, somewhat lower 
than the model predicted value of 0.5±0.2. In four cases, 



one can measure also the ratio of the 11.3 to the 6.2 /mi 
feature. The values are about 0.5 in 3 cases, and 0.8 (with 
an unceratinty of ±0.2) in the forth object, in good agree- 
ment with the template model predictions of 0.4±0.2. It 
is likely, as we will discuss further in §5.3, that our models 
tend to overestimate the intensity of the 3.3 /im feature 
by a factor 2-3 in stars with 3 x 10 4 < x < 3 x 10 5 . 

However, there are two groups of non-detections that 
are interesting. We discuss these ones in the following. 

5.1.2. Flat disks 

The first is the group of 5 diamonds, i. e., HAeBe stars wit h 
no evidence of flared disk (group II. lMeeus et~aTl l200l|) . 
Of these, only in 1 case (HD 142666) PAHs have been 
detected at 3.3 and 6.2 /tm; most of the non-detections 
are well below the model predictions. This agrees quite 
well with our result that flat disks should have very weak 
PAH features, even when PAHs are present with normal 
properties on the disk surface. 

5.1.3. High-x stars 

Other interesting objects are those with high values of 
X (> 3 x 10 5 ). Most of them (the exception being WL 
16, t hat we discuss in mo re details in the following) are 
from IBrooke et all (^93), and only the 3.3 /mi feature 
has been observed. All these objects with high x have 
3.3 /im intensities (or upper limits) well below the model 
predictions. 

The interpretation of this group of objects is not 
straightforward; it is quite possible that their disks, if they 
ever existed, hav e been dissipated by the strong radiatio n 
field of the star llNatta et all 1200(1 iFuente et all 12002^ . 
or that they are flat l|Hillenbrand et allll992h . There are, 
however, other possibilites. One is that the abundance of 
PAHs is lower than we assumed. Photoevaporation, which 
increases rapidly with x and for smaller PAHs, could eas- 
ily reduce the PAH abundance well below our assumed 
value. PAHs evaporation, as discussed in A3. 31 is a com- 
plex process, and we may have underestimated its effects. 

It is also possible that the PAHs are ionized and/or de- 
hydrogenated in the inner disk region, resulting in a much 
reduced 3.3 /tm band (see t!4.2.2l and l4.2.1JI . Observations 
of the PAH emission from WL 16, shown by a square 
in Fig. OH 4 , favors this last possibility. In fact, we find 
that while the models reproduce well the 11.3 /tm fea- 
ture strength, that of the 3.3 /mi feature is well below the 
predictions. Moreover, IRessler fc Barsonvl l)2003|) interpret 
the spatial variation of the 7.7 (or 8.6)/11.3 /tm flux ratio 
as due to a change of the PAH charge state along the disk. 
In the inner region PAHs appear to be positively charged 

4 In Fig. [5J the data for the embedded young stellar object 
WL 16 has been corrected f or dust attenuation taking a visual 
extinction of Av = 31 from Ressj^^^^arson^ 12003!) and us- 
ing th e extinction coefficient of iDram^^^Leg (^84fi ; iDrainel 
(1989). 
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while in the outer region they are neutral. Furthermore, 
the 12.7/11.3 ^m flux ratio which increases in the outer 
part of the disk could indicate the presence of more hy- 
drogenated PAHs in the outer part. 

5.2. Which PAHs? 

As we have just seen, our template model, with relatively 
large (Nc = 100) neutral PAHs fits very well the strength 
of the most commonly observed features, with the possible 
exception of the 3.3 fim one, which tends to be much lower 
than predicted in stars with x>3x 10 5 . 

In stars with lower x (spectral type roughly later than 
B9), we can rule out that PAHs are large and ionized, be- 
cause they will produce too low 3.3/6.2 and 11.3/6.2 /im 
ratios (by a factor of >10 and 4, respectively, see H4.2.1I 
and TableOJ). Similarly, large and strongly dehydrogenated 
PAHs are unlikely, since they will also produce too low 
3.3/6.2 and 11.3/6.2 ^m ratios (see H4.2.2JI . Further, we 
can rule out that PAHs are much smaller than our as- 
sumed model (No = 100), because they will produce too 
high 3.3/6.2 and 11.3/6.2 /im ratios (by a factor of 4 and 
10, respectively, see ij4.2.3l and Table 0). 

On the other hand, we cannot exclude that PAHs are 
small and ionized. In this case, we find that the strength 
of the 6.2 and 11.3 /jm features, as well as their ratio 
will be roughly similar to those predicted by the template 
model with large and neutral PAHs (see t!4.2.3J) . This is 
easily understood considering that the effects of the size 
and ionization compensate each other. For small positively 
charged PAH, the compensation is nevertheless not com- 
plete, and in fact, the intensity of the 3.3 fim feature will 
be lower than in the template model by a factor ~10 (see 
Table EJ). Therefore, this model, which underestimates the 
3.3 fim feature strength by a factor of 3 in stars with 
3 x 10 4 < X < 3 x 10 5 but by a much larger factor in 
stars with x < 3 x 10 4 , seems unlikely. For small nega- 
tively charged PAHs, the 3.3 /jm feature is expected to be 
stronger (see i!4.2.1|) and could become comparable to the 
observations. However, it is difficult to assess the reliabil- 
ity of this model, since the IR properties of PAH anions 
are not well known. 

To get some additional insight into the PAH prop- 
erties, one can look at the other weaker PAH features 
(i.e., 8.6, 12.7 or 16.4 fim) even if they have been only 
sporadically detected. This is probably due to the low 
S/N ratio of the spectra. The models with large neutral 
PAHs and small ionized PAHs predict similar strengths 
of the 8.6 and 12.7 fim features, in good agreement with 
observations (see Tables and |2J). On the other hand, 
large neutral PAHs will produce a 16.4 /im feature five 
times stronger than small ionized PAHs. This favors the 
small ionized PAHs model since this feature has not 
been detected in the ISO spec tra of isolated HAeBe stars 
ijVan Kerckhoven et all l2002|) . However, high S/N spec- 
tra in this wavelength region are needed to determine the 
properties of the carrier of this band. 



The large and neutral PAHs hypothesis is supported 
by the detailed studies of the PAH spectra in HAeBe 
stars base d on both laboratory data and theoretical cal- 
culations jffonv et allEioTllVai! Kerckhoven et alll2002t 
IPeeters et all 12002^1 . Together with the comparison of 
the PAH profiles observed in evolved stars and ISM-like 
sources, they indicate that PAHs should be larger and 
less ionized in isolated HAeBe stars, in agreement with 
our template model. 

It must be emphasized that there are several compli- 
cations which we have neglected in this study. The most 
obvious is that we have assumed that PAHs can be char- 
acterized by a single size and charge state. This is unlikely 
to be the case, and one can expect variations as a func- 
tion of radius and depth in the disk, as well as from object 
to object. For example, PAHs are likley to be more pos- 
itively ionized in the inner disk region and for stars with 
higher \- This could explain the gradual decrease of the 
3.3 /im feature strength with x seen in Fig. El In stars with 
3 x 10 4 < x < 3 x 10 5 , it is possible that we have an inter- 
mediate situation, and, in fact, a model with a mixture of 
both large neutral and small ionized PAHs will predict a 
3.3 /im feature strength lower than in the template model 
by a factor of 2-3, close to the observed values. 

We conclude that in stars with x < 10 5 PAHs may be 
preferentially relatively large (Nc >100) and neutral. But 
small and negatively ionized P AHs can also reprod uce the 
observations. For comparison, iLi fc Luninel l)2003h . which 
have modelled the emission features from the (debris) disk 
around HD 141596A (with x = 10 5 ) conclude that PAHs 
must be mostly negatively charged. However, they show 
that models consisting of a mixture of both neutral and 
charged PAHs are also capable of reproducing the obser- 
vations. For stars with x > 10 5 , the large neutral PAH 
model tends to overestimate the intensity of the 3.3 fim 
feature. This could result from changes in the disk struc- 
ture (e.g., disks may have dissipated), or from changes in 
the characteristics of PAHs, which can be ionized, dehy- 
drogenated or strongly photoevaporated in the inner disk 
regions. 

From the integrated spectra, it is difficult to discrim- 
inate between the various effects due to ionization, de- 
hydrogenation or evaporation. In order to provide bet- 
ter insight into the PAH properties evolution within the 
disk, high angular spatial spectroscopic observations in the 
bands are needed. 

6. Summary 

We have investigated the emission from PAHs on the sur- 
face of disks around HAeBe stars by comparing model pre- 
dictions with observations. We have computed models of 
disks, heated by irradiation from the central source, which 
contain large grains, in thermal equilibrium with the ra- 
diation field, and transiently heated very small grains and 
PAHs. The disks are optically thick to the stellar radia- 
tion, in hydrostatic equilibrium in the vertical direction, 
with dust and gas well mixed (flared disks). We have used 
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a 2-layers model to calculate the disk structure and imple- 
ment a ID radiative transfer code to compute the emerg- 
ing spectrum. This scheme is reliable and efficient (see §3.1 
and 3.2). Our main results can be summarized as follows. 

1. The models predict an infrared SED showing PAH fea- 
tures at 3.3, 6.2, 7.7, and 11.3 /zm clearly visible above 
the continuum. The PAH emission, spatially extended, 
comes mostly from the outer disk region (R ^100 AU) 
while the continuum emission at the same wavelength, 
mostly due to warm large grains, is confined to the in- 
nermost disk regions (R ~few AU). Among the PAH 
features, the 3.3 /im one is the least extended; in our 
template model, about 1/2 of its integrated intensity 
comes from R <30 AU. 

2. From comparison with ISO and ground-based obser- 
vations, we find that most of the obser ved PAH fea- 
tures in objects with flared disks (group I. lMeeus et all 

are well described by our mo dels. In objects with 
no evidence of flared disk (group II. lMeeus et allEoOlh 
no (or weak) PAH features have been detected, as ex- 
pected from our models. For geometrically flat disks, 
PAH features are predicted to be very weak, even when 
PAHs are present with standard properties on the disk 
surface. 

3. Objects with strong radiation field (spectral type gen- 
erally earlier than about B9) have the 3.3 /jm feature 
(often the only one observed) much weaker than pre- 
dicted. We suggest that their disks have been dissi- 
pated by the strong radiation field, or that PAHs have 
been photoevaporated more effectively than our model 
predicts, or that their properties (ionization or dehy- 
drogenation) have been modified in the inner disk re- 
gion. All these scenarios will cause the 3.3 /im feature 
to practically diseappear. 

4. Finally, we find that in circumstellar disks around stars 
with spectral type later than B9, PAHs should have 
an abundance typical of the ISM. They may be pref- 
erentially large and neutral, or small and ionized, and 
integrated spectra alone cannot distinguish between 
these two possibilities. In order to provide better in- 
sight into the PAH properties evolution within the 
disk, high angular spatial spectroscopic observations 
are needed. The fact that we can explain the observed 
PAH spectra, at least of the isolated HAe stars, with 
PAH abundances and qualitative properties similar to 
those of PAHs in the ISM is interisting. It seems that, if 
PAHs are reformed in disks, after having been depleted 
(by coagulation or other processes) in dense cores, they 
contain again roughly the same amount of carbon (rel- 
ative to that of large grains) that they contain in the 
ISM. This is a very interesting point, which certainly 
deserves further attention. 

To the best of our knowledge, this study represents 
the first successful modeling of the PAH emission features 
from disks around pre-main-sequence stars. It shows that 
PAHs are present in disks, at least around HAeBe stars, 



and absorb a significant fraction of the stellar radiation, 
similarly to PAHs in the ISM. Consequently, PAHs are 
an important source of opacity in circumstellar disks and 
are likely to play a dominant role in the thermal budget 
and chemistry of the gas. The coupling between PAHs and 
gas can determine the gas temperature in the outer disk 
regions, where the lower density can make the collisional 
coupling between gas and grains inefficient, and as a con- 
sequence have a significant impact on the disk structure 
itself. Finally, as a side aspect, PAH emission which is sen- 
sitive to the local physical conditions - either directly or 
through a change of their physical characteristics - could 
be useful to better understand the physical conditions pre- 
vailing in the surface of circumstellar disks. Furthermore, 
PAH emission, which is more extended than the thermal 
emission at similar wavelengths, can offer an unique pos- 
sibility to probe the outer disk structure. 

Acknowledgements. We are grateful to F. Boulanger and L. 
Verstraete for fruitful discussions and their relevant comments 
and suggestions. 
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Appendix: Numerical method for the radiation trans- 
fer 



The disk is divided into rings of radius r and small 
radial width Ar. At the center, at r = 0, sits the star 
with parameters M*, L*, R* and T e ff. For each ring, we 
compute, as described below, the radiative transfer per- 
pendicular to the disk, in z-direction, and then add up 
the contributions from all rings. 

Each ring is treated as a plane-parallel slab. It displays 
mirror symmetry with respect to the mid plane at z = 0. 
We solve the integral equations 



I+(t) = e- T (V(0) + jT S{x)e x 



dx 



J-(r) = e" T Jo + 



S{x) e x dx 



(1) 

(2) 



I + (r) is the intensity at optical depth r and frequency 
v in upward direction along a straight line under an angle 
p = cos 9 with respect to the z-axis; r is zero at z = 0. The 
second equation containing I~ refers to downward beams. 
Iq is the intensity of radiation that falls from outside (the 
star) on the disk. The indices 9 and v have been dropped 
for convenience of writing. S(t) is the source function. The 
optical thickness r and geometrical height z are related 
through dr = p k dz where p is the mass density and k 
the mass extinction coefficient of dust. 

There are two boundary conditions. One states that 
at the top of the slab the incoming intensity I~ = I is 
zero except towards the star where it equals -B„(T c ff). The 
stellar hemisphere, which seen from a distance r subtends 
a solid angle £1 = 7ri?J/2r 2 , is replaced by a luminous 
band of the same intensity B v (T e g). The band encircles 
the sky at an elevation a and has a width Q/2n <C a, its 
solid angle is also O. The other boundary condition reads 
1+ = I- at z = 0. 

One starts the calculation of the radiative transfer for a 
ring with a guess for the vertical dust temperature profile 
T(z). This fixes the source function S(t) once the dust 
properties have been specified. Then one solves Q and 
J5J for I+(z,p) and I~(z,p) and calculates from I + ,I~ 
the mean intensity of the radiation field J(z). This allows 
to compute a new temperature run T(z). The procedure 
is iterated until T(z) converges. 

Because in our examples the visual opacity of the disk 
in z-direction can be very large (> 10 5 ), we apply the fol- 
lowing slight modification to what we said above. We slice 
the disk into a completely opaque mid layer sandwiched 
between two thin top layers. The mid layer, because of its 
high optical depth and the lack of internal energy sources, 
is considered to be isothermal at temperature Tmid! the 
value depends, of course, on r. We evaluate the radiative 



transfer for the top layer according to the method de- 
scribed before, and separately for the isothermal midlayer, 
but there it is trivial. A boundary condition connects the 
two layers. In the calculation for the top layer, we now 
position the zero point of the z-axis at the transition be- 
tween top and mid layer. The intensity that enters the top 
layer from below under the direction p = cos 9 is 



I+{z = Q,p) = B„{T mid )(l-e- T ^ 



(3) 



Here t„(/i) is the optical depth of the mid layer. The 
monochromatic flux received by the top layer from below 
is therefore 



F v = 2-kB v {T mid ) J (l-e- T »^pdp 



(4) 



The frequency integral J /„ dv must be equal to the 
downward flux at z — because the net flux, for reasons of 
disk symmetry, vanishes everywhere. This condition yields 
T m id, of course, iteratively. In mathematical form, it reads 



dvB u (T mid ) j dfi(l-e- T "^ 



dv 



dp pl v (z = 0) (5) 



Ideally, the top layer extends downwards so deep that 
the temperature at its bottom asymptotically approaches 
the temperature of the mid layer, T m i<j. Then the whole 
disk of arbitrary thickness is calculated consistently. But 
that would require a visual optical depth of the top layer in 
z-direction, T t0 p, of several hundred and therefore many 
(hundreds of) vertical grid points. In practice, we have 
found that T top ~ 4 is sufficient, since pratically all 
starlight is absorbed and T(z — 0) is very closed to T m i d . 
Typically, changing r t0 p from 4 to 40 changes the emerging 
SED by less than 20% at all wavelengths. 

We wish to make a few further points: a) The radia- 
tion that is absorbed in the top layer and then reemitted 
at infrared wavelengths may leave the top layer, be again 
scattered or absorbed there, or may enter the mid layer. 
b) The net flux is always zero, but the upward and 
downward integrated fluxes, F + and F~ , do depend on z, 
so always F + + F ~ =0, but dF + /dz ^ 0. c) Good agree- 
ment in the results of our radiative transfer code with 
another one is demonstrated under the internet address 



http: / /www. mpa-garching.mpg.de/dullemon/radtrans/benchmarks 
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Fig. 9. Comparison of model predictions and observations. Upper Panel : Strength of the 3.3 /im feature normalised 
to the FUV flux of the star for ZAMS stars. As in Fig. the scale on the upper axys is the spectral type, that on 
the lower axys is the value of x, *-e-, the FUV field at 150 AU from the star. The two solid lines show the predictions 
of the template model for <fi+ = 0.5 (lower curve) and (/)+ — 1.0 (upper curve) (see text). The two thin lines show 
an uncertainty strip of ± a factor of 2 around these models. The crosses show observed intensities (normalized to 
the corresponding FUV flux) for the isolated HAeBe stars with evidence of flared disk, the diamonds for the isolated 
HAeBe stars with no evidence of flared disk and the triangles fo r the intermediate s ources (HD 97048 and Elias 1). 
The asterisks show the HAeBe stars observed from the ground bv lBrooke et all IWO.i ): the square is WL 16, corrected 
for dust attenuation (see text). Error bars, when not shown, are smaller than the symbols. Arrows are 5a upper limits. 
Medium and lower Panels : Same as upper pannel for the 6.2 and 11.3 fim feature. 
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Table 1. Parameters of Herbig Ae/Be candidates. Evidence for a flared disk and disk outer radius from spectral energy 
distribution modeling and millimeter interferometry observations. Presence of AIBs and silicate bands. Crystalline 
signatures. 



Object 


d 


Sp. Type 


T eff 


L* 


a 

X 


Flared 


r> b 


AIBs 


Sil. 


Cryst. 


Ref. 




[pc\ 




r ts~\ 


T 

L© 




disk 


[AU\ 










AB Aur 


144 


B9/A0Ve 


9750 


47 


10 b 


J 

V 


140 


J 

V 


>/ 

V 




(1) 


HD 100546 


103 


B9Ve 


11000 


36 


10 5 


V 


300(150 bl ) 


v 7 


V 




(1) 


HD 142527 


200 


F7IIIe 


6250 


31 


4 10 3 


V 




yj 


V 




(1) 


HD 179218 


240 


B9e 


10000 


80 


2 10 5 


V 


(<600 62 ) 


V 


V 




(1) 


HD 100453 


114 


A9Ve 


7500 


9 


5 10 


V 




V 


~ 


" 


(1) 


TTTi 1 O f O A A 

HD 135344 


84 


F4Ve 


6750 


3 


6 10 2 


/ 

V 










/ l \ 

(!) 


HD 139614 


157 


A7Ve 


8000 


12 


10 


/ 

V 










(!) 




1 A K 
140 


/YO ve 


1UOUU 




i n ,r > 

1U 


/ 

V 


zzu 


/ 
V 






/ 1 "i 


HD 104237 


116 


A4Ve 


10500 


40 


10 5 


110 






J 

V 


? 


(1) 


HD 142666 


116 


A8Ve 


8500 


11 


10 4 


110 




V 


V 


•v/ c 
V 


CD 


HD 144432 


200 


A9Ve 


8000 


32 


3 10 4 


no 






V 


? 


(i) 


HD 150193 


150 


AlVe 


10000 


40 


10 5 


no 






J 

V 


V 


(i) 


HD 163296 


122 


A3Ve 


10500 


30 


10 5 


no 


360 




V 


V 


I'D 


HD 97048 


180 


B9-A0 


10000 


31 


10 5 




f30o 6i -Qoo b2N i 


V 






(2) 




150 


A0-A6 


8000 


21 


2 10 4 




70 ( <r?nn b 3 \ 


,/ 

V 


,/ 

V 


? d 


(2) 


T t Hrw 9^ 
1 J K MM ■ ) 




An 

nil 


youu 


Q 


9 1 n 4 




f^l 900"! 
^ V* izuu ^ 


/ 
V 






{<>) 


I 1 t\ WI 1 




A r i 


OJUU 


1 9 


in 4 

1U 






./ 

V 






{<>) 


tin Ozl^l 8^ 

III.' i'iJIOO 


Ann 


A 9 


qi nn 

yiuu 


1 7 

1 1 


^ i n 4 

O 1U 




~\n\ 
\ lw ) 


,/ 

V 






{<>) 


T k Hn 1 Q8 


Ron 


A 


ooou 


R 

u 


in 3 




(Qnn\ 




,/ 

V 




(■'<] 

v-v 


BD+61°154 


650 


B8 


11200 


330 


10 6 




(900) 








(3) 


T Ori 


460 


B9 


10700 


83 


3 10 5 




(88) 








(3) 


V380 Ori 


460 


B9 


10700 


85 


3 10 5 




(150) 








(3) 


HD 250550 


700 


B7 


12300 


240 


10 6 




(410) 








(3) 


Lk Ha 208 


1000 


B7 


12300 


210 


10 6 




(1600) 








(3) 


VV Ser 


440 


B9 


10700 


100 


4 10 s 




(410) 








(3) 


WW Vul 


700 


A3 


8700 


20 


3 10 4 




(88) 








(3) 


BD+46°3471 


1000 


AO 


9800 


508 


10 6 




(190) 








(3) 


Lk Ha 233 


880 


A5 


8300 


30 


3 10 4 




(320) 








(3) 


WL 16 


125 


B8-A7 


9000 


240 


4 10 b 




(450) 








(4) 



^ : detection, ?: possible detection, -: no detection. 

a Far-ultraviolet (FUV, 6< hu <13.6 eV) flux at 150 AU from the star expressed in units of the average interstellar radiation 
field, 1.6 10" 6 W in" 2 faabind Il968h . 

In parenthese we giv e esti mate of the PAH emission extension. References: (bl) Ivan Boekel et al] §004); (b2) 
Siebenmoreen et alJ (|200CT) : fb3l iBrooke et alJ lll993ft . 
c Possible blend with the PAH 11.2(3) /jm band. 
d lHanner et"aT1 dl994Tl. 

References: (D lMeeus etal] J200lf) and refere nces therein; (2~) IVan Kerckho ven^^J] (|2002|) and references therein; ( 3) for astro- 
physic al parameters iHille'nbrand^et alJ (Il992l) (except for BD+46°3471 see iBerrillietalJ al992h and for VV Ser seejTesti^^al 



1998)); fo r PAH emission extension a nd AIBs detection lBrooke et alJ il993h : for silicate emission features detection lBerrilli et al 
1992 ); f4) lRessler fc Barsonvl i2003T) and references therein. 
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Table 2. Integrated strength of the PAH emission features. 



Object 


3.3 


6.2 


"7.7" 


8.6 


11.2(3) 


12.7 


16.4 


Ref. 




H 


\fj,m] 


[H 


H 










AB Aur 


<i a 


4.4(0.3) 


4.2(0.6) 


2(0.2) 


2(0.6) 






(1) 


HD 100546 


2.5(0.5) 


14.3(0.4) 


19.2(0.8) 


4.6(0.5) 








(1) 


HD 142527 


1(0.6) 


<1.5 






<0.5 






(1) 


HD 179218 


1.7(0.2) 


8.1(0.4) 


14.4(0.6) 


1.6(0.4) 




- 


- 


(1) 


HD 100453 


1.3(0.2) 


2.6(0.5) 


~ 


_ 


1.2(0.3) 


- 




(1) 


HD 135344 


<0.5 


<1.5 


? 




<1 


- 




(2) 


HD 139614 


<0.5 


<1.5 


~ 




<0.5 






(2) 


TTT\ 1 /iA-| a r\ 

HD 169142 


1(0.2) 


2.8(0.3) 






2.3(0.5) 






(!) 


HD 104237 


<0.5 


<1.5 












(2) 


HD 142666 


0.3(0.2) 


2.6(0.3) 












(1) 


HD 144432 


<0.5 


<1.5 












(2) 


HD 150193 


< 0.08 a 


<1.5 












(2) 


HD 163296 


<0.2 a 


<1.5 












(2) 


HD 97048 

111.' {j 1 Ut:U 


1.3(0.3) 


6 Q(T) 8"! 


11.8(1.3) 


2.4(0.6) 


3.4(0.4) 


1 1'D 2) 






J_JlldO _L 


n c i c fn D 




Q 4(0 1 ^ 
ij .'-tyxj. _L ^ 




d 






("PI 


T k Hrv 9^ 


n n7fi(Ti nnfil 

U.U I U^U.UUU J 
















tri/ r\ r : 
nrv w 1 1 


U.Uo / ^U.U^iO ) 














a\ 
[■') 


nn 94^1 sec 


















ijiv net iyo 


•rn i ^ 
















DLJ-\-\) 1. lO'rt 


vu.uy 














{■>) 


T Ori 
















\°) 


V380 Ori 


<0.1 
















HD 250550 


<0.12 














(3) 


Lk Ha 208 


<0.03 














(3) 


VV Ser 


<0.05 














(3) 


WW Vul 


<0.09 














(3) 


BD+46°3471 


<0.1 














(3) 


Lk Ha 233 


<0.04 














(3) 


WL 16 


0.17(0.03) 








5(0.4) 






(4) 



Integrated fluxes (after continuum substraction) and uncertainty (in between brackets) are in units o f 10~ 14 W/m 2 . For non- 
detections, we give for the 3.3, 6.2 and 11.3 urn feature 5cr limi ts estimated from the spectrum of iMeeus et al] ll200l]) and 
assuming constant feature widths taken from lLi fc Drainel (l2001al) . We have not included the 11.2(3) band for the sources 
with st rong silicate emission . 
a From lBrooke et all (Il993h . 

6 In the spectrum of HD 1005 46 a 11.2(3) /jm PAH feature is present but may be blended with the crystalline olivine band at 
11.3 /an iMalfait et al.lll998l) . 

c Measured in ground -based spectrum llGeballeLll997h . 

d lHanner et al] ( Il994f) report the presence of an 11.2 /jm feature in their UKIRT spectrum. They attribute this feature to the 
presence of crys talline silicate but it could als o be a PAH band. 

References: (1) IVan Kerckhoven et al.l <2002ft: (21 IMeeus et all ]200ll) ; (3) iBrooke et all ]l993Tl : (4) From the spectrum of 
iTokunaea et al 1 <199lfl and lHanner et al] il992^ for the 3.3 and 11.3 fim feature, respectively. 

Table 3. Predicted strength of the PAH emission features integrated over the disk. 



Model 


3.3 


6.2 


7.7 


8.6 


11.3 


12.7 


16.4 




H 


H 


H 


[urn] 


H 


H 


H 


template 


2.8 


5.4 


10.2 


1.3 


2.5 


0.6 


2.6 


Ionized 


0.1 


7.5 


18.2 


2.25 


0.9 


0.13 


0.8 


/h=0.5 


1.4 


5.9 


11 


1 


1.5 


0.3 


2.6 


A c =40 


5.3 


2.5 


4.1 


1.1 


10.8 


2.8 


1 


JV C =40 - Ionized 


0.3 


6.1 


12 


3 


3.9 


0.9 


0.4 



Integrated fluxes (after continuum substraction) are in units of 10~ 14 W/m 2 . 



